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Cluster structure of O, O and O is investigated by the antisymmettrized molecular dynamics 
(AMD) plus generator coordinate method (GCM). We have found the K'^ and Oj" rotational 
bands of that have the prominent ^''C+a cluster structure. Clustering systematics becomes 
richer in ^°0. We suggest the K^=0'^ band that is the mixture of the ^^C+a+4n and ^^C+^He 
cluster structures, and the K^ =Q~[ band that has the '^''C+^He cluster structure. The K^—Q^ and 
0^ bands that have the prominent ^^C+o? cluster structure are also found. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

The investigation of the cluster structure of atomic 
nuclei has been mainly performed for light N=Z nu- 
clei, and much less is known for N^Z and unstable 
nuclei. A famous example of the clustering in N^Z 
nuclei is the 2a clustering of neutron-rich Be isotopes 
[H, 0, i, i, i, S 0, i . The AMD studies 0, i, 0, Q suc- 
cessfully described many properties of neutron-rich Be 
isotopes and showed that most of them have the 2a clus- 
ter core. The degree of 2 a clustering dynamically changes 
depending on the motion of valence neutrons. This fea- 
ture is well understood by the concept of the molecular- 
orbital P, 0, Q. The AMD study % has predicted the 
existence of the cluster structures also in ^■^Ne, that is 
a system expected to have the ^^0-\-a cluster core as an 
analogue of neutron-rich Be isotopes. In this study, it has 
been showed that the ^^0-|-a core is formed and dissolved 
depending on the configurations of two valence neutrons. 
Thus, the cluster structures of these nuclei suggests us 
richer variety of the cluster structure in A^t^Z nuclei than 
N=Z one. 

To understand the clustering in N^Z nuclei, a series 
of Oxygen isotopes is another good start point, because 
the "'^^C-l-a cluster structure in -'^^O has been investigated 
in detail for a long time [lol |. Our interest in this study 
is to understand what will happen to the ^^C-|-a cluster 
structure when we add neutrons to ^^O. 

Many of the low- lying excited states of ^^0 are un- 
derstood within the shell model space of two neutrons 
in sd-orbital above N—Z—8 shell closure. However, it is 
also known that the states with the core excitation (4p2/i) 
coexist [U in the same energy region. They were as- 
sociated to the ^"^C-f a cluster structure, and many works 



were performed to investigate the molecule band struc- 
ture of [H d [ll d [H 113, [S H [131. Conse- 
quently, the K'^=0~^ band that consists of the observed 
0+ (3.63 MeV), 2+ (5.26 MeV), 4+ (7.11 MeV) and 6+ 
(11.69 MeV) states has been established well as the 4p2/i 
and ^"^C-l-a molecular band. The cluster model calcula- 
tions [lB[3 tl^3,t showed the ^"'C+a cluster structure of 
the K'^=0^ band also predicted two molecular K^=0~ 
bands. More recently, the a-cluster structure was in- 
vestigated by the elastic a scattering on [ll] and 
a breakup reaction [H, [H, [H]. In Ref. [23|, it was 
proposed that the K^=0~ band built on the 1^ state 
at 8.04 MeV [2^, that was in same energy region with 
the K'^=02 band predicted by the cluster models, was 
the parity doublet partner of the a-cluster band. 

In the case of ^°0, nothing is known about the cluster- 
ing. In the investigation of the neutron transfer reac- 
tion ^^0(t,p)^°0 compared with the shell model predic- 
tions, it has been suggested that the core excited states 
(6p2/i) coexist with the normal states (ApOh) We 
consider that these states (0+(4.46 MeV), 2+ (5.30 MeV), 
and 4+ (7.75 MeV)) could be associated with the cluster 
state as in the case of ^^O. 

The purpose of this study is to investigate the clus- 
ter structure of ^^O and ^"O. We have applied the 
AMD-I-GCM (antisymmetrized molecular dynamics plus 
generator coordinate method) framework. The AMD is 
a kind of ab initio theory in the sense that it can describe 
the cluster structure and shell-like structure within the 
same framework without such assumptions as existence 
of clusters or an inert core. Therefore this framework is 
useful to investigate the existence of the cluster struc- 
ture in ^^O and ^"O. We also calculate ^^O to investigate 
how the ^^C-|-a cluster states are described in this frame- 



2 



work. We will suggest various kinds of cluster states in 
and 20O. We find the K^'^O'^ bands that have promi- 
nent ^"^C+q; cluster structure, and some negative parity 
states that have non-negligible amount of ^^C-l-a com- 
ponent. In ^"^O, the motion of valence neutrons around 
the ^^C+a cluster core enriches the variety of the clus- 
tering. Depending on the motion of valence neutrons, 
-^^C-|-a-|-4n, ^'^C+^He and ^^C+a cluster structures ap- 
pear. 

The contents of this article are as follows. In the next 
section, the AMD-I-GCM framework is briefly outlined. 
In Sec. mil the cluster structures of ^^O, and are 
discussed. In the last section, we summarize this work. 



II. THEORETICAL FRAMEWORK 

A. AMD wave function and calculational 
procedure 

In this subsection, the AMD-I-GCM framework is 
briefly outlined. For more detail, readers are directed 
to Refs. [13, m . The AMD intrinsic wave function of 
A-nucleon system is described by a Slater determinant. 
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(1) 
(2) 



Here, fi is a single particle wave packet which is com- 
posed of spatial part (pi (r) , spin part Xi i and isospin part 
Ti. The spatial part is described by a Gaussian, 
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(3) 



where is complex three dimensional vector. The width 
parameter ly is common for all nucleons and fixed to 0.17 
fm^^. Spin part is parameterized by complex number 
parameter ^i, 



(4) 



The isospin part is fixed to up (proton) or down (neu- 
tron). The Zi and are the variational parameters and 
optimized by the frictional cooling method. The parity 
projected wave function which is generated from the $int 
is used as variational wave function. 



where the Px is the parity operator. 

The Hamiltonian used in this study is given as. 



H = T + Vn + Vc-T, 



(5) 



(6) 



The f is the total kinetic energy and Tg is the energy 
of the center-of-mass motion, that is exactly treated in 



the AMD. As the effective nuclear force Vn, the Modi- 
fied Volkov force (MVl) [1^ and spin-orbital part of the 
G3RS [3^ force are used. Details and the applied param- 
eter set of these forces are given in the next subsection. 
Coulomb force Vc is approximated by a sum of seven 
Gaussians. 

The energy variation is performed under the constraint 
on the matter quadrupole deformation parameter (3. The 
constraint potential. 



Vcnst = VcnstiiP) - Pof ^ 



(7) 



is added to the total energy of the system. Here, Vcnst 
takes adequate positive value, and /3o is a given number. 
The definition of {(3) is given in Ref . [3| • 

After the variation, the optimized wave function (/3) 
is projected to an eigenstate of the total angular momen- 
tum J, 

The integrals over three eular angles are calculated nu- 
merically. 

Finally, we superpose <I'^^(/3) and diagonalize the 
Hamiltonian. The wave function which describes a cer- 
tain state is given as. 



(9) 



where Ci is determined by the Hill- Wheeler equation, 

5{{K^\H\K^) - 6„(*,i±|*„'±)) = 0. (10) 

B. Interactions 

We use the MVl caseS force [ll] for central force, and 
the G3RS force ^ for spin-orbit force. The MVl force 
consists of finite-range two-body and zero-range three- 
body terms. 



=1/(2) +f(3)^ 

^ - m - mP„Pr) 



(11) 



x{V:4exph [h./TAf] + VrgM- (ryAfl)']},(12) 



i<j<k 



(13) 
(14) 



The spin-orbit part of the G3RS force is given as, 
Vls = J2^{^"''''' - e~""^-}P('0)iy • (s, + Sj),(15) 



where P('^O) is the projection operator onto the triplet 
odd state. Adopted force parameters are summarized 
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in Table. U and the binding energies of ^He, i244(-i 
and 16,18,200 calculated using the present AMD+GCM 
framework are shown in Table, [ill For ^^C, the excitation 
energies of the 2f and 4^ states which are considered to 
be important to describe the low- lying states of ^^O are 
also shown. 



TABLE I: The force parameters of the MVl case3 force and 
spin-orbit part of the GSRS force. 



m Va [MeV] 


Vr [MeV] rA [fm] 


rR [fm] [MeV] 


0.61 -83.34 


104.86 1.60 


0.82 4000 


u [MeV] 


Ki [fm^^] 


tin [im^''] 


3000 


5.0 


2.778 



TABLE II: The binding energies of i6,i8,20q^ 4jjg 12^ 
The excitation energies of 2]*" and 4^ states of ^^C are also 
shown. 





B.E [MeV] (Cal.) 


B.E [MeV] (Exp.) 


*He 




28.9 


28.29 


"C 




102.8 


105.28 


160 




127.3 


127.62 


ISq 




137.5 


139.81 






153.2 


151.36 


@ 




B.E [MeV](0+) 


Ex. (2+) [MeV] Ex. (4+) [MeV] 




Cal. 


89.4 


5.5 13.3 




Exp. 


92.16 


4.44 14.1 



C. Analysis of the single- particle orbits 

We also investigate the single-particle structure of 
the obtained wave function <I>^(/3) by diagonalizing the 
single-particle Hamiltonian Q. First, we transform the 
single-particle wave packets ipi into the orthonormal basis 

'fa, 

1 ^ 

iPa ^ —=^Cia(pi. (16) 

Here fia and Cia are the set of eigenvalues and eigenvec- 
tors of the overlap matrix Bij = {Lpi\ifj), 

A 
i=i 

Using the (pa, we construct the single-particle Hamilto- 
nian matrix, 

A 

haf3 = {(pa\i\^p) + '^^{(pa(f>i\v\(pi3'f>t - tfi^Pp) 
i=l 

1 



Then we obtain the single-particle energy ep and single- 
particle wave function by the diagonalization of h^fj, 

hapgiip = tpQap, (19) 

P 

A 

0s = ^ gap<fa- (20) 

In this study, we calculate the density distribution of the 
single-particle wave function (f>s to investigate the motion 
of valence neutrons. 



III. RESULTS 
A. I'^O 

is well known to have the prominent ^^C-f a cluster 
structure in its excited states [1^. First, we investigate 
1^0 and see how the cluster structure is described in the 
present framework. 

The energy curves before and after the angular mo- 
mentum projection for the (a) positive- and (b) negative- 
parity states are shown in FIG. [T] Before the angular 
momentum projection, the positive-parity curve (dotted 
line in FIGlT](a)) has a energy minimum at the spherical 
point. As the deformation becomes larger, the energy 
rapidly increases. The angular momentum projection 
drastically changes the energy curve. Here, we discuss 
each curve with respect to the angular momentum eigen 
states with K—0 for the sake of simplicity. The O'^ curve 
has a minimum at /3=0.20 and a shallow local minimum 
at /3=0.66. The minimum state has the Ohuj configura- 
tion and correspond to the ground state, while the local 
minimum state has the AHuj configuration (proton 2huj 
and neutron 2huj) and contributes to the state. Here 
the particle-hole configuration of each state is evaluated 
by the analysis of single particle orbits. The density dis- 
tributions of the intrinsic wave functions at these minima 
are shown in FIG.[2](a) and (b). As clearly seen, the wave 
function at /3=0.66 has the prominent ^^C+a clustering. 
In the case of the 2+, 4+ and 6+ curves, they have two en- 
ergy minima around /3=0.30 and 0.65. The energy curve 
of the negative-parity state is also steep before the angu- 
lar momentum projection. After the angular momentum 
projection, the energy spectra are different in the moder- 
ately deformed region {f3 < 0.5) and the largely deformed 
region {(3 > 0.5). In the former region, the lowest state 
is the 3~ state, and 1~ state is approximately 5 MeV 
above the 3~ state. In the largely deformed region, the 
spectrum shows the rotational nature. It is due to the 
structure change of the intrinsic wave function. In the 
moderately deformed region, the wave function has the 
Ihuj configuration, while in the largely deformed region, 
it has the ^^C+a cluster structure as shown in FIG.[2](c) 
and (d). 

After the angular momentum projection, we have per- 
formed the GCM calculation. The states with non-zero 
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K quantum number are also included into the GCM cal- 
culation. FIG. [3] shows the calculated and observed level 
scheme. We have obtained the excited K'"=GI and 
rotational bands together with the ground state and the 
low-lying 3j~, 1^ and 2^ states. The ground state dom- 
inantly consists of the wave functions around the mini- 
mum at /3=0.20. The excited K^'^Q^ and Oj" rotational 
bands mainly consist of the wave functions in the largely 
deformed region (/3=0.5-0.7 and /3=0.6-0.8, respectively) 
that have the prominent ^^C-|-a cluster structure. The 
low-lying 3^, 1^ and 2^ states consist of the wave func- 
tion that have the Ihw configuration. The excitation en- 
ergy of the 3j~, 1^ and 2^ states and the moment of 
inertia of the excited K'"=Q'l and bands are quali- 
tatively reproduced. However, the calculated energy of 
the K'^=0'l and 0^ bands considerably overestimate the 
experimental value. We consider that one of the rea- 
son is that the internal wave function of ^^C(g.s) cluster 
is not correctly described in the present framework. In 
our wave function that has the ^^C-fa cluster structure 
(FIGHKb) and (d)), the ^^C is not in its ground state but 
has some excited state component. This difficulty would 
be overcome if we performed the double angular momen- 
tum projection (a simultaneous projection with respect 
to the ^■^C cluster and to the total system), though it 
is rather time consuming. To compare with the exper- 
iment, we utilize the ^^C-|-a threshold energy evaluated 
without the projections. Namely, we evaluate it from the 
energies of '^He and ^^C calculated without the parity, an- 
gular momentum projection and GCM, and call it theo- 
retical threshold energy in the following. The theoretical 
^^G+a threshold energy is 17.1 MeV (the energies ^^C 
and ''He are —83.5 and —26.7 MeV, respectively), and 
it is shown in FIGl3]by the dotted line. When their en- 
ergies are measured from the threshold, the K'^=Q'^ and 
Oj~ bands show better agreement with the experiment. In 
the present results, the parity doublet bands constructed 
by the prominent ^^C-t-a cluster structure are consistent 
with the cluster model calculations In addition, the 
single-particle excitations in the low-lying negative parity 
states are also described. 



In this subsection, we investigate how the a cluster 
structure changes by adding two neutrons to -'^^O. The 
same calculational procedure as the case of ^^O is ap- 
plied to ^*0. The obtained energy surface and the den- 
sity distributions of the core and valence neutrons are 
shown in FIG. [4] and [5l respectively. Here, we have de- 
fined the valence neutrons as two neutrons in the most 
weakly bound neutron orbitals, and the core as the nucle- 
ons in the lowest 16 orbitals. In the states shown in FIG. 
[5l two valence neutrons occupy the orbitals that have the 
same spatial density distributions. The 0^ curve has the 
energy minimum at /3=0.20, and the 2+ and 4"*" curves 
have energy minima in the /3 ~ region. These minimum 
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FIG. 1: Energy curves of ^''O as functions of matter 
quadrupole deformation parameter /3 for the (a) positive- and 
(b) negative-parity states. Solid line represents the energy of 
each parity and angular momentum (if^O) state, and dashed 
lines show the energy before the angular momentum projec- 
tion. 
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FIG. 2: The matter density distributions of the intrinsic wave 
functions of ^^0. The centroids of the single-particle wave 
packets are plotted with white squares. The "l'i„t(+) and 
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denote the intrinsic wave function on the positive and 



negative parity curves, respectively. 
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FIG. 3: The excitation energies of the low-lying states of 
^^0. Energies of the excited sates are shifted, for ease of com- 
parison. Dotted lines show the theoretical and experimental 
threshold energies. 



states have the Ohuj configuration, although the density 
distribution of the 0"'' minimum state (FIG. [5] (a)) shows 
small deformation with the parity asymmetry. Around 
/?=0.45, the 6"*" curve has the energy minimum, and the 
0"*", 2+ and 4+ curves have the shoulder. In this re- 
gion, the wave functions are approximately correspond 
to the proton 2Haj configuration. The density distribu- 
tion FIG. [5] (b) shows that the system is separated into 
two clusters. There are 14 wave packets in the left side 
and 4 in the right side, that indicates the formation of the 
^"^G+a cluster structure. Indeed, the density distribution 
of two valence neutrons shows that the valence neutrons 
stay only around the ^^C cluster. The formation of the 
^^C+a cluster structure leads to the rotational nature 
of the 0+, 2+, 4+ and 6+ energies. The wave functions 
around /3=0.51 become the dominant component of the 
K^^02 rotational band after the GCM calculation. 

In the case of the negative parity states (FIG. ID^b)), 
the 3~ curve has the energy minimum at (3—0.23, where 
the wave function has the proton Itiuj configuration. The 
1~ curve has the energy minimum at /3=0.34. The den- 
sity distribution of this state (FIG.[5](c)) shows the shght 
development of the cluster structure. As deformation 
becomes larger, this cluster structure develops. Figure 
[5] (d) shows the pronounced ^"'C+a cluster structure of 
the largely deformed negative parity state. Although the 
energy curves have no local minimum, they become the 
dominant component of the K^=Q^ rotational band after 
the GCM. Again, the development of the ^'^C-\-a clus- 
ter structure is confirmed from the distributions of the 
wave packets and the localization of the valence neutrons 
around ^^C (FIG. [5] (d)). 

The level scheme obtained by the GCM calculation is 
shown in FIG. [Altogether with the experiment. We have 
obtained the ground band and many excited states in- 
cluding the K'^=02 and 0^ bands that have the "'^''C-l-a 
cluster structure. The ground band (the 0^, 2^ and A'^ 
states) dominantly consists of the wave functions around 
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FIG. 4: Energy curves of ^**0 as functions of matter 
quadrupole deformation parameter /? for the (a) positive- and 
(b) negative- parity states. Notations are same with FIG. [T] 
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FIG. 5: The density distributions of the core (black contour 
lines) and the valence neutrons (color plots) of ^*0. The cen- 
ters of the single-particle wave packets are shown by the white 
squares. The intrinsic wave function of (a) and (c) gives the 
minimum energy for the 0^ and 1^ state, (b) and (d) be- 
comes the dominant component of the /S''^=Oj and K^=Q'^ 
rotational bands. 
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/3=0.20. Compared to the experimental spectrum, the 
level spacing in the ground band is considerably under- 
estimated. It is mainly due to the strong spin-orbit in- 
teraction used in this calculation. When we weaken its 
strength as u=2000 MeV, the 2^ and 4^ states come to 
better positions {2f state at 1.0 MeV and 4:^ state at 2.0 
MeV). However, the use of the weaker spin-orbit interac- 
tion leads to the overestimation of the ^'^C+a threshold 
energy. Therefore we use the present strength parameter 
to discuss the ^''C-l-a structure. 

We have found that the Oj, 2;^, 4 J and 6^ states dom- 
inantly consist of the wave functions around /3=0.51, and 
hence classified them as the K'^—O'^ band. As mentioned 
above, the intrinsic wave functions around (3—0.51 have 
the ^'^C+a structure, therefore this band is regarded to 
have the ^'^C+a cluster structure. The ^"^C+a threshold 
energy is overestimated as in the case of ^^O. Therefore, 
to compare the obtained i^'^=Oj rotational band with 
the experiment, we measure its excitation energy from 
the threshold energy. The theoretical ^''C-l-a threshold 
energy calculated in the same way as the ^'^C+a thresh- 
old energy is 10.3 MeV (the energy of is —100.5 
MeV), while the experimental one is 6.24 MeV. The en- 
ergies of the 0^, 2'^, ^rid 6^ states measured from the 
"C-Ha threshold energy are -2.1, -0.6, 0.7 and 4.8 MeV, 
respectively. They approximately agree with the exper- 
imental 0+ (3.63 MeV), 2+ (5.24 MeV), 4+ (7.11 MeV) 
and 6+ (11.69 MeV) states that are at -2.61, -1.00, 
0.87 and 5.45 MeV measured from the threshold energy, 
respectively. This assignment is consistent with that pro- 
posed in many theoretical works [111, H, [iM and the 
a-transfer reaction j2l|. We stress that the ^'^C+a clus- 
ter structure of the if =0^ band has been found without 
assuming any structure in our calculation. We note that 
there is the mixing between the ^^C+a cluster structure 
and the shell-like structure. The if^=2^ band (the 2^, 
3^ and 4^ states) mainly consists of the wave functions 
around /3=0.20, that have the Ohuj configuration. The 
natural parity states with this configuration are mixed 
with the ^'^C+a cluster wave functions {K=0). There- 
fore the 2^ and 4;^ states also have considerable amount 
of the ^"^C-l-a component, although the amount of their 
cluster component is smaller than that of the ii"^=Oj 
band. 

The K'^=0^ rotational band (Ig , 3^ and 52^) domi- 
nantly consists of the wave functions around P—0.72, that 
have the prominent ^"^C-l-a cluster structure. Therefore, 
the can be regarded as the parity doublet part- 

ner of the K'^~02 band. The energies of the 1;^, 3^ and 
5^ states measured from the theoretical ^'^C-l-a threshold 
energy are 2.5, 3.0, and 7.28 MeV, respectively. Experi- 
mentally, the assignment of the K'^—0~ band which con- 
sists of the 1- (8.04 MeV), 3- (9.67 MeV), and 5" (11.62 
MeV) sates was proposed by the a breakup reaction . 
The energies of these states measured from the ^''C-l-a 
threshold energy are 1.80, 3.43, and 5.38 MeV, respec- 
tively. We consider that they are candidates of the cal- 
culated K'^=Oi band. As in the case of the positive par- 




FIG. 6: The low-lying level scheme of ^ O. The experimental 
candidates of the K^=02 and 0^ bands are quoted from [2ll . 

m 

ity states, there is the mixing between the ^'^C+a struc- 
ture and the proton Ifko configuration. This leads to the 
fragmentation of the ^'^C+a cluster structure into many 
states. Especially the lj~, I2 , 3j^, 5^ and I4 states have 
considerable ^^C+a cluster state component, though it 
is much smaller than that in the K'^=0~ band members. 



By adding four neutrons to ^^O, a variety of cluster 
states appears. The obtained energy curves and the den- 
sity distributions of the core and four valence neutrons 
are shown in FIG. [7] and [U respectively. We define the 
valence neutrons as the four neutrons in the most weakly 
bound neutron orbitals, and the core as the nucleons in 
the lowest 16 orbitals. In the states shown in FIG. [HI 
there are always two orbitals that have different density 
distributions, and two valence neutrons occupy each or- 
bital. The 0^ curve has the energy minimum at P—0.20 
(FIG. [8] (a)), and the 2+ and 4+ curves have the energy 
minimum around P—0.10 and /?=0, respectively. In this 
region, the wave functions have the OSw configuration. 
Around P—0.41, the structure changes from the Ohuj to 
the proton 2hu! configuration, and the O"*" and 2+ (4+ and 
6^) curves have the shoulder (local minimum). In this 
region, two different cluster structures appear. Let us 
compare the wave functions at (3=0.41 and (3=0.53 (FIG. 
[8] (b) and (c)). Both of them have similar core density 
distributions which show the development of the ^^G+a 
cluster core. The difference between them is clearly seen 
in the density distributions of the valence neutrons. The 
density distribution FIG. [8] (b) shows that four valence 
neutrons orbit around entire ^'^C+a core. On the con- 
trary, FIG. [H] (c) shows that two of four valence neutrons 
localize around ^^C cluster and the others localize around 
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FIG. 7: Energy curves of ^"O as functions of matter 
quadrupole deformation parameter /3 for the (a) positive- and 
(b) negative- parity states. Notations are same with FIG. [T] 



a cluster. Therefore, we regard that the wave function 
at P—OAl has the ^^C-fa-|-4n structure, and the wave 
function at (3—0.53 has the ^^C-f^He structure. In the 
former, valence neutrons are moving in the mean field of 
the whole system of the ^^C+a core, while in the latter, 
the spatial correlations of two neutrons with the ^^C and 
a core are enhanced. From /3=0.58, where the 0+, 2+, 
4+ and 6+ curves have shoulder, another cluster structure 
appears. The density distribution of the wave function 
at /3=0.62 (FIG.[8](d)) shows the formation of developed 
^^C-|-a structure, in which all valence neutrons orbit only 
around the ^^C cluster. 

Various structures also appear on the negative parity 
curve. The 1~, 3~ and 5~ curves have energy minimum 
around /3=0.17, where the intrinsic wave functions have 
the proton Ifko configuration. Around (3—0.50, where 
the negative parity curves show the rotational nature, 
the ^''C+^He structure appears. The density distribu- 
tions of the core and valence neutrons in this state (FIG. 
|5] (e)) are quite similar to those of the ^'^C-l-^IIe struc- 
ture that appears on the positive parity curve (FIG. [H] 
(c)). In the largely deformed region, the ^^C+a cluster 
structure appears around (3=0.61 (FIG. [5] (f)), which is 
quite similar to that found on the positive parity curve 
(FIG.Kd)). 

The low- lying level scheme of obtained by the 



GCM is shown in FIG. [9] together with the experiment. 
The ground band (K'^^Ot), many excited rotational 
bands {K''=Oj, Oj, 2'^ , 2^, 0^ and 0^) and other ex- 
cited states have been obtained. The level spacing in 
the ground band (0^, 2^ and 4^ states) is considerably 
underestimated compared to the experiment. It is due 
to the strong spin-orbit interaction, as in the case of 
1^0. The K'^=0^ band (Oj, 2+, 4+ and 6^ states) dom- 
inantly consists of the wave functions around /3=0.41- 
0.58. In this region the ^^C-|-a-|-4n and ^''C-l-'^He clus- 
ter structures appear, as mentioned above. Therefore, 
the K'^=02 band is the mixture of these structures. If 
we consider that the calculated state correspond to 
the experimental (4.46 MeV) state, the experimen- 
tal 2+ (5.23 MeV) or 2+ (5.30 MeV) state and 4+ (7.75 
MeV) states are the candidate of the calculated /C'^=Oj 
band from their energy positions. In the experimental 
side [2I, these 0+ (4.46 MeV), 2+ (5.30 MeV) and 4+ 
(7.75 MeV) states have been assigned to the proton 2huj 
states predicted by the analysis with shell model calcula- 



tions. The K'^=22 band also consists of the wave func- 



The K''=0^ band (the I3 , 3^ 



4 states) dominantly consists of the wave functions 



tions around /?=0.41 
and 5^ 

around (3=0.50, that have the ^''C-l-^He cluster structure. 
Since the ^'^C-l-^He cluster state component is contained 
also in the K'^=02 band, the K^=02 and Oj~ bands can 
be interpreted as the parity doublet bands, although the 
K^=02 band has some ^^C-|-Q;-f4n cluster state compo- 
nent. In the negative parity states, there is the mixing 
between the proton Ihoj and the ^^C-|-^He state as in 
the case of ^*0. This results in the fragmentation of the 
i4(-;_|_6jjg cluster state component into the I4 , 34 and 
53" states. Therefore these states have non-negligible 
amount of the ^"^C-l-^He cluster state component. The 
K^=0^ band dominantly consists of the wave functions 
around /3=0.62, that have the ^^C+a cluster structure. 
The K^=02 band also consists of the wave functions 
that have the ^^C+a cluster structure. As already men- 
tioned, the ^^C-I-q; structures in the positive and neg- 
ative parity states around /3=0.62 are quite similar to 
each other, therefore these bands are regarded to be the 
parity doublet bands. These bands start from just above 
the theoretical ^^C+a threshold energy. However, the ex- 
perimental information on the levels around the ^^G+a 
threshold (12.32 MeV) is very little unfortunately. The 
K^=2^ band also consists of the wave functions around 
(3=0.62. 

Let us discuss characteristics of the cluster features of 
^"O in a series of O isotopes. As mentioned above, we 
have found that the a-cluster states follow the threshold 
rule systematically in the series of ^^O, ^^0 and ^"O. 
Namely, the K'^=0^ band of ^^O {^^C+a structure), 
K-^=0+ band of ^^O {^^C+a structure) and K''=0^ 
band of (^^C-f-a structure) appear near the corre- 
sponding threshold energy. On the contrary, the K'^=02 
band of ^''O appears at much smaller excitation energy 
than the "C-H^He threshold energy, although this band 
has large component of the ^''C-f^IIe cluster structure. 
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FIG. 8: The density distributions of the core and the valence neutrons of O. Two valence neutron orbitals are shown for each 
intrinsic state. 



We note that the ^-^C-l-^He wave function (FIG.[5](c)) has 
large overlap with the ^^G+a+An wave function (FIG. [5] 
(b)), and hence the system has molecular-orbital-like na- 
ture. Therefore, we consider that the valence neutrons 
play an important role to lower the energy of the K^=02 
band. The presence of the molecular-orbital-like band in 
^"O may be related to the weakly bound nature of ^^C 
and ®He. Since the last two neutrons in both nuclei are 
weakly bound, ^^C-I-q; and ^"'C+^Hc do not appear in 
small excitation energy. However, when two neutrons 
are covalently bound and shared by ^^C and a clusters, 
it lowers the energy of the system. 



IV. SUIVIIVIARY 

We have investigated the cluster structures of ^^O, ^^O 
and •^"O, using the AMD-I-GCM framework. First, we 
have confirmed that the K^=Of and Oj~ bands of ^^O 
have the ^^C+a cluster structure. In ^^O, The K^'^O^ 
and 0^ bands that have the ^^G+a cluster structure are 
obtained. They appear around ^'^G+a threshold energy 
and can be regarded as the parity doublet bands. It is 
also noted that the ^''C-l-a cluster structure is fragmented 
into many states. The valence neutrons give richer struc- 
ture for ^°0. The analysis of the valence neutron or- 
bitals revealed the presence of the cluster structures that 
have different motion of valence neutrons. First is the 



"'^^C-|-Q;-|-4n cluster structure, in which four valence neu- 
trons orbit entire ^^G+a core. Second is the ^''C-l-^IIe 
cluster structure, in which the valence neutrons localize 
either of the ^^C or a cluster. These structures are mixed 
and construct the K^=Q2 band. This band has much 
smaller excitation energy than the ^"^C-f^He threshold 
energy. The ^"^C-l-^He cluster structure also constructs 
the K^=Q^ band that can be regarded as the parity 
doublet partner of the K'"=Q^ band. Third is the ^^G+a 
cluster structure that constructs the K'^=Q^ and 0^ par- 
ity doublet bands around the ^^C-fa threshold energy. 
The appearance of the variety of cluster states in 
may be related to the weak binding nature of subsys- 
tems, ^^C and ^He. They are not rigid cluster subunit, 
because the last two neutrons are weakly bound in both 
nuclei. As a consequence, four valence neutrons orbit 
around ^^C-f a cluster core in different ways, and it leads 
to the variety of cluster structures. It is an open prob- 
lem whether the clustering phenomena appears in further 
neutron-rich Oxygen isotopes. 
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